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Abstract 

Based on the AdS/CFT correspondence, string theory has given exact predictions for circular 
Wilson loops in U(N) Af = 4 supersymmetric Yang-Mills theory to all orders in a 1/N expansion. 
These Wilson loops can also be derived from Random Matrix Theory. In this paper we show that 
the result is generically insensitive to details of the Random Matrix Theory potential. We also 
compute all higher fc-point correlation functions, which are needed for the evaluation of Wilson 
loops in arbitrary irreducible representations of U(N). 
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1 Introduction 



One of the surprising predictions of the Maldacena conjecture (T[ relating type IIB string theory on 
an AdSz, x S 5 background to M = 4 supersymmetric Yang-Mills (SYM) theory in four dimensions 
concerns the static quark-antiquark potential. At strong 't Hooft gauge coupling A = g 2 N this 
potential is predicted to be of strength \f~X~ [||, |3| rather than A (as at weak coupling). Because the 
gauge coupling does not run in M = 4 SYM this could in principle be tested by explicit gauge theory 
computations. Substantial progress in this direction was made last year by Erickson, Semenoff and 
Zarembo Q . They computed the sum of all Feynman-gauge planar diagrams without internal vertices 
to the expectation value of a circular Wilson loop in N = 4 SYM and found (I n {x) is the nth order 
modified Bessel function) 

(W) circle = ^=h(Vx) ~ sJl^i , (l) 

in agreement with the leading-order (in A) AdS / CFT prediction |5|, [6| 

(W) AdS/C FT = . (2) 

By an explicit computation it was also shown in ref. 0] that the first corrections to ([!]) of order A 2 
coming from diagrams with internal vertices precisely cancel in four dimensions. 

Recently, Drukker and Gross Jj], |8| have made a quite remarkable extension of this result by (a) 
outlining a proof that the above sum of rainbow diagrams (|l|) actually gives the exact result to all 
orders in y/X at N = do, and (b) showing that the calculation can be extended to all orders in a 1/N 2 
expansion. Comparing again with the AdS/CFT correspondence they find exact agreement to leading 
order in A, at every order in 1/N 2 . 

One of the key ideas of ref. |7j is to make efficient use of a Random Matrix Theory (RMT) 
interpretation pj] of the result (Q). Consider the unitary ensemble of N x N hermitian matrices, and 
a corresponding partition function Z of Gaussian Boltzmann factor. Identifying N with that of the 
gauge group U (N) of J\f = 4 SYM, it was noticed in H that in the limit N —* oo 



— TV e M exp 

N 



A 



= ^i(VA), (3) 

precisely reproduces the gauge theory result (||). Intuitive reasons for why this is a correct represen- 
tation of the circular Wilson loop in M = 4 SYM have been given in refs. [f|, 0. Moreover, in ref. 
Q it was argued that this particular Random Matrix Theory representation can be used to compute 
exactly the leading order (in vA) circular Wilson loop to all orders in the expansion. It should 

thus be possible to completely by-pass the otherwise quite cumbersome computation in the N = 4 
supersymmetric gauge theory language. 

The fact that the RMT representation was chosen with Gaussian Boltzmann weight is clearly 
linked directly to the fact that only rainbow graphs are being summed in the field theory context. An 
obvious question to ask is why this can yield the full answer. Although arguments have been given 
0, f7j, a full proof is still lacking. While the resulting effective theory of the circular Wilson loop can 
be argued to be zero-dimensional (basically on account of a "conformal anomaly" invalidating the 
conformal mapping from the line to the circle at just one space-time point 0), there is a priori no 
guarantee that vertices cannot contribute at that particular point. Such interaction vertices would 
correspond to higher order terms in a more general RMT potential V{M). 

The purpose of the present letter is two- fold. First, we shall turn the question around, and ask for 
the consequences of including higher-order terms in the RMT potential V{M). Remarkably, there is 



2 



a huge degree of universality at work here. To leading order in 1/N 2 the precise result (||) changes, 
but the leading- order contribution in \f\ remains generically unaffected^ Moreover, all higher-order 
terms in 1/N 2 are generically insensitive to the precise form of the potential V(M), but depend only 
on a finite number of "moments" M k . These statements of universality follow from a series of results 
derived in connection with the RMT approach to 2D quantum gravity || [HJ . Our second purpose is 
to advocate the use of the so-called loop insertion method [11] to compute these circular Wilson loop 
expectation values. We believe this technique is superior to the orthogonal polynomial method in the 
present context, independently of whether one wants to restrict oneself to the Gaussian potential or 
not. To illustrate this, we show how to compute the most general A;-loop correlation function for any 
potential V(M), to all orders in the 1/N 2 expansion. We first briefly, in the section below, recall some 
of the main results from Random Matrix Theory. 

2 The /c-point function to all orders in 1/N 2 

Our starting point is the observation that the loop expectation value ((l/iV)Tr exp[y/\M]) is the 
Laplace transform C of the 1-point function or resolvent G(p). Defining 



ZT'ICIriKx) = (^Tre 1 ") = W(x) . ('.) 



we immediately have 

rrMcUnWf^ = / 

N 

Here, the expectation value with respect to the Random Matrix Theory partition function of the 
unitary ensemble 

r 00 l 

Z = dMexp[-NTrV(M)} , V(M) = ^ v9kM k . (6) 
J k=i k 

is defined with a general potential V(M) in the standard way. 

To compute a Wilson loop in an arbitrary representation of U(N) we need also more general 
expectation values 

±Tv e XlM ■ ■ — Tr e x " M ) . (7) 



,N N 

These can be computed from the corresponding connected loop expectation values 

00 -. 

W( Xl ,...,x k ) = N k - 2 (Tre^ M ---Tre^ M ) = ^ Wfa, . . . , x k ) . (8) 

\ / conn. £ — ' JS J 

j=0 

The latter will follow similarly from the connected /c-point correlation function (/c-point resolvent) 

1 _ 1 \ A 1 

conn. j—Q 



G(p u ..., Vk ) = N k - 2 (Ti -•••Tr -) = ^ ^(p!,...,^) (9) 



by k inverse Laplace transforms. For that purpose we will rely on the very powerful results for the 
connected /c-point function which have been given explicitly in ref. [11] to all orders in the 1/N 2 



(genus) expansion. The reader familiar with this material can skip it, and go straight to section 3. 
1 We shall comment on the precise condition below. 
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2.1 The 1 -point function 

Starting with the 1-point function the general solution is of the following form [11] 



3A;— 1 

gh(p) = E (4V°(p)+^V w (p)) 

n=l 



fc > 1. 



(10) 



The constants A^ 1 and D^ 1 are rational functions of the "moments" 



M k = 
Jk = 



do; 



duj V'juj) 

h%iri(w-ay/ 2 (Lj-b) k +V 2 



9k+i + 9k+2 Q& + {k + ^)aj + ■■■ 



9k+i + 9k+2 [ -a + 



+ 5 »] + ... 



ii 



and depend on at most 2(3fc — 1) of them. These moments encode in a universal way the dependence 
on the infinite set of coupling constants {gk}, and we will give a few examples below. The functions 
X (p) an d ip( n \p) are the basis functions needed for the solution of the loop equation. They are given 
explicitly in ref. plj[ . All we need to know here is that x^ n \p) an d ip^ n \p) are linear combinations of 
the functions 4>a{p) = {p— a)^ k ^ 1 / 2 (p — b)^ l l 2 and (fy^ (p) = (p— &) -fe-1 / 2 (p — a)" 1 / 2 , respectively, up 
to order k = n. They also depend on the moments and up to k = n. Since the inverse Laplace 
transforms of </>a^(p) an d <\>l (p) ar e known in closed form, we immediately read off the corresponding 
expansion of the Wilson loop W(x). For a Gaussian potential V(M) = \g2M 2 we have b = —a and 
Mi = J\ = 4/a 2 = gi- All higher moments M/%>2 and Jfc>2 vanish. In particular the correlators then 
contain less terms and thus differ from the general expression. 

The result given in eq. (|lC|) gives only the 1/N 2k corrections to the 1-point function G(p) for 
k > 1. The leading order result of k = is non-universal and depends explicitly on all the coupling 
constants in the potential V(M). It can be written |l2| 



G (p) 



1 



du> V'{uj) (lu — b)(u> — a) 



+ 



1 



2 Jc 2m p — u y (p — 6) (p — a) ^/(p ~~ — a ) 
The boundary conditions that fix the end-points of the cut [b, a] are 



hk = -{ ^ ; 1 



fc = 0,1 



(12) 



(13) 



which follow from the requirement that Go(p) ~ | for p — > oo. In order to illustrate the non- 
universality let us give an explicit example for a symmetric sixth-order potential: 



G (p) 



1 



92P + 54P 3 + 9eP 5 - ( 92 + ^a 2 + ^9(>a 4 + (54 + ^g&a 2 )p 2 + c/ 6 p 4 ) \jp 2 - a 2 



The endpoint a of the cut (which is now symmetric) [—a, a] is given by eq. (13): 



— floQ H Q4CL H Ofia . 

4 y 16 y 32 yb 



(14) 



(15) 



2.2 Higher k-point functions 



Higher /c-point functions can be obtained by applying successively the loop insertion operator [11] 

d ^ k d 



dV(p) 



E 

k=l 



p k+1 dg k 



(16) 
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starting from the resolvent itself: 

G(pi, 



d 



d 



,Pk) 



dVipu) dV{p 2 ) 



G(pi) • 



(17) 



The planar 2-point function obtained in this way from eq. ( |12D is universal, the well-known result of 
Ambj0rn, Jurkiewicz and Makeenko 0]: 



G (p,q) 



A{p-qf 



(p - b)(q - o) + (p - a)(q - b) 
V(P -b)(p- a)(q -b)(q- a) 



(18) 



As a consequence all fc-point functions for k > 2 are universal to all orders in 1/N 2 , because the 
derivatives d/dV(p) acting on the universal variables b,a,Mk,Jk can be expressed again in terms of 



these variables and the functions 4^ip) an d <fi'a > {p)- They are again given explicitly in |Q]. 

Apart from Go(p,q) all universal Appoint functions Gj(p±, . . . ,pk) factorize with respect to their 
arguments into a finite number of terms to all orders in 1/N 2 . This follows from the fact that once 
a function Gj(p%, . . . ,Pk-i) is factorized, the application of the loop insertion operator d/dV(pk) only 
adds factors which are functions of p&. At genus zero the 3-point function factorizes 



(*)/ 



G {p,q,r) 



1 



1 



1 



1 



Ml (p-a) 3 / 2 (p-6)V2 (g_ a )3/2( g _ 6)1/2 ( r _ a )3/2( r _ 6)1/2 
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Jl (p-a)V2(p_ 6)3/2 ( g _ a )l/2( g _ 6 )3/2 ( r _ a )l/2( r _ 6 )3/2 



• (19) 



By applying successively the loop insertion operator one sees that all Go(pi, . . . ,Pk) also factorize for 
k > 3. At higher genus factorization follows immediately from the 1-point function eq. (10) (see 
the example for genus one below). This factorization property makes the successive application of k 
inverse Laplace transformations particularly simple. The only exception is Go(p,q), which requires a 
little more work. 

Let us finally mention that the same procedure we have described here applies immediately to 
the complex matrix model |L3| ] and to the so-called reduced hermitian matrix model. There, a closed 
expression has been derived for all planar Appoint correlators with k > 2 || 14]. All higher genus 
corrections in 1/N 2 are also available for both the complex [15] and reduced hermitian models [Hi 



Precisely as for the hermitian case treated below, we can simply read off the corresponding loop 
correlation functions from the inverse Laplace transformations. 



3 Wilson loops from Laplace transforms 

We are now ready to compute Wilson loop expectation values by means of inverse Laplace transforms. 
We begin with the Wilson loop in the fundamental representation of U(N). As already mentioned, this 
expectation value as computed in Random Matrix Theory is not universal. Consider as an example 
the sixth-order (symmetric) potential from eq. (|i~4|). To leading order in 1/N 2 we get 

2 3 / 5 \ a 3 5 a 5 

W (x) = —h(ax) + -\gi + -g&a 2 — I 3 (ax) + — g 6 — I 5 (ax) . (20) 
ax 8 V 4 / x 32 x 

This shows that the natural variable is the combination ax, and we identify a = The asymptotic 
behavior for large A appears to depend in a complicated way on the coupling constants g\ and gg. But 
from the constraint (15) it follows that in order to achieve a — > oo we require at leastf^s^n ~ a~ 2n (or 



2 For an arbitrary symmetric potential the constraint ( |l~3| ) can be written 1 = 5 2~ZtLo 9 2k Cfc ) TS" 
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even higher suppression in 1/a) . Because the leading term in the asymptotic Bessel- function expansion 
at large argument is I n (z) ~ exp [z] / V 2irz for all fixed n, the large-A behavior is therefore 

W (x) ~ const. A~ 3/4 exp[x\/A] , (21) 

in agreement with the Gaussian result [Q, fjj when we set x = 1. So the leading large-A behavior is 
unaffected by higher order terms in the potential V(M). 

We now turn to the first 1/N 2 correction to this result. For simplicity, let us here again restrict 
ourselves to symmetric potentials (but see the Appendix for the most general expressions). The 
universal function W\ follows from the genus-one 1-point function G\(p) for a symmetric potential 

r ( x 1 aM 2 1 , , 

XVP) 4Mi (p 2 - a 2 ) 5 / 2 8(Mi) 2 (p 2 - a 2 ) 3 / 2 ' 1 J 

Taking the inverse Laplace transform, we get the general result 

W 1 (x) = ^L- * 2 / 2 M - ^ x h (ax) . (23) 

Specifically, for the same symmetric sixth-order potential as above we have 

3 15 

Mi = 92 + ^a 2 g 4 + — a 4 g 6 , M 2 = 2ag A + 5a 3 g 6 . (24) 
z o 

For a Gaussian potential this reproduces the result of ref. (7|: VKi(A) = A/2(v / A)/48. Again we find 
that the leading large-A behavior is unaffected by higher order terms in the potential, after using 
M\ ~ a~ 2 ,M2 ~ a~ 3 as a — > oo: W r i(A) ~ const. X 3 ^ exp[\/A]. Note that the second term in (p3|) is 
subleading in this limit. We have also computed the general genus- two (1/N ) contribution G2(p), 
and confirmed the result of ref. 0] for the special case of a Gaussian potential. Clearly, the leading 
A-behavior is also here of the same form for any generic potential. 

We next turn to higher fc-point correlation functions of these fundamental Wilson loops, needed 
for the evaluation of Wilson loops in arbitrary representations of U(N). The universal 2-loop result 
(|i"8|) is seen to not factorize in p and q. One may use the convolution theorem to evaluate the needed 
double inverse Laplace transform, but it is simpler to note that 

d 1 pq ~t~ oP 

Q^G (p,q) = 4 (p2 _ a2)3/2((?2 _ a2) 3 /2 (25) 

does factorize. It is thus elementary to get the inverse Laplace transforms, 

—W (x,y) = ^[I (ax)I (ay) + h{ax)h{ay)} (26) 

which we integrate up to 

W Q (x,y) = ^ f du u[Io(ux)I {uy) + h(ux)h{uy)] , (27) 
2 Jo 

after fixing the integration constant by comparing [|| 

G °(PiP) = = Za IT + ° ( Zz) = 4(TrMTrM) conn . + O ( ±) (28) 



4(p 2 -a 2 ) 2 p 4 4 \p 6 J p 4N ' \p' 

and 

W (x,y) = xy (TrMTrM) conn . + O (x 2 y 2 ,xy 3 ,x 3 y) . (29) 
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Performing the integral in eq. ( |27|) we get for an arbitrary symmetric potential 

W (x,y) = QXy [I (ay)h(ax) + I (ax)h(ay)} . (30) 
2(x + y) 

As a check on this result, we note that it vanishes when one or both of the arguments x and y become 
equal to zero. Finally, at equal arguments we have 

ax 

W (x,x) = —I (ax)I 1 (ax) . (31) 
As an example, let us evaluate A r_2 ((TVexp[M]) 2 ), which in ref. (7| is denoted by W± : i. We get 

= (W (l)) 2 + ^[W (lA)+2W (l)W 1 (l)}+O^-^j , (32) 
where all functions on the right hand side have been given above. In particular, we find 



in the Gaussian case. The leading large-A behavior of the 1/N 2 correction goes like exp[2vA], which 
is now ready to be compared with string theory through the AdS/CFT correspondence. 

It is straightforward to go on to arbitrarily high order in 1/N 2 for any n-point correlation function. 
As a last example, consider the universal 3-point function Wo(p,q,r) for an arbitrary symmetric 
potential. We find: 

W (x,y,z) = xyz^-[I 1 (ax)I (ay)I (az) + I (ax)I 1 (ay)I (az) 

+I (ax)I (ay)Ii(az) + I 1 (ax)I 1 (ay)I 1 (az)\ (34) 

where for the Gaussian case M\ = gi = 4/a 2 = 4/ A. The most general result for the genus expansion 
of the resolvent can be written 

3fc-l 

W k (x) = (4" ) ^ 1 [X ( ' ,) (P)]W+4 T,) ^ 1 ^ ( ' ,) (P)]W) . *>1- (35) 

n=l 

where the inverse Laplace transforms are known explicitly. Any n-point function for n > 2 can then 
be found to all orders in the l/N 2 expansion by the iteration procedure described above. We have 
thus succeeded in determining all n-point Wilson loop correlation functions to all orders in the 1/N 2 
expansion. They are all given by appropriate products of n modified Bessel functions. 

4 Other universality classes 

By fine-tuning the RMT potential one can reach new (multi-critical) universality classes that would 
invalidate the above conclusions. To see what happens in such cases, it suffices to focus on the 1-point 
function, the Wilson loop itself, in the planar limit. First note that, with p(Q indicating the eigenvalue 
density 



^Tre xM \ = 1 JdM ^Tr e xM exp [-NTrV(M)) = £ d& exp(a<i) (36) 



where for simplicity we have restricted ourselves to even potentials V(M), and the cut has been 
rescaled to lie on the interval [-1,1]. By inserting the Wigner semi- circle law p{Q = (2/tt)(1 - Q 2 ) 1 ' 2 
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corresponding to a Gaussian potential we of course just recover the result © with a = VX (The 
form ( |36| ) also immediately shows the non-universality of this result, away from the large- A limit). 
Multicritical densities belonging to multicritical universality classes near the soft edge can be chosen 

Pm(o ~ a - c 2 r +1/2 ■ (37) 

According to eq. (36) the corresponding mth multicritical Wilson loop is 



(k + lV2 k+1 e xy ^ 
Wo{x) = (ax) k +i 4+l(ax) ~ COnsL A (2fc+3)/4 ' ( 38 ) 

The leading large-A behavior thus differs in the prefactor from that of generic non-critical potentials, 
although the exponential form ~ exp[x\/A] is retained. We have explicitly checked eq. ( |38|) for m = 1, 2 
by tuning the potential accordingly in eq. (|14|). For this special choice of coupling constants the Bessel 



functions in eq. ( |20|) precisely combine to give eq. (|3q). 

It is tempting to speculate that the universality classes of multicritical points could correspond 
to circular Wilson loops of different gauge theories at conformal points. But we have no evidence to 
support such a claim. 
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A Appendix 



In this appendix we give, for completeness, a few examples to illustrate that the method we have 
presented also works for an arbitrary potential that is not restricted to be even. We start with a 
closed expression for the non-universal Wo(p) which can be derived from eq. (|l~2| ) 



W (x) 



+ exp x 



du exp u 



a + b 



a + b 



If) u 



In X 



a — b 



(lo — b)(uj — a) 
(p ~ b)(p - a) 



exp[u>(x — u)] 



(39) 



Next, we find the 1/N 2 contribution: 



Wi{x) 



16 



exp x 



16 



exp x 



) \M 1 
1 

+ n 

a + b\ |Y M 2 



\{M 1 



12 1 

-I {xc) + -h (xc) + -h (xc) 

lib 

12 1 

-J (xc) - -h (xc) + -I 2 (xc) 

' X [J (xc) + h (xc)} 



+ 



cM x 



(40) 



where we have defined c= (a — b)/2. Note the appearance of the exponential prefactor exp[x(a + 6)/2] 
which disappears for the symmetric potential with b = —a. The result for a generic symmetric 
potential eq. ( |2"3| ) is easily recovered by setting = (— l) k+1 Mk and c = a. As a last example let us 
also give the planar connected 3-loop function 



W (x,y,z) 



a — b 
xyz — - — exp l xyz- 



^r) \tt t J o (xc) + A (xc)] [J (yc) + h (yc)] [J 



I Mi 



[zc) + Ji (zc) 
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j- [J (xc) - h (xc)} [I (yc) - h (yc)} [I (zc) - h (zc)} }. 

(41) 



Also the general expression for the 2-loop function can be evaluated for general asymmetric potentials, 
using a variant of the method described above for the symmetric case. Proceeding iteratively, this 
again determines all n-point functions to all orders in the l/iV 2 -expansion, now for potentials that are 
not necessarily symmetric. 
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